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Abstract 

In this paper we investigate tours of the Jovian satellites Europa, Ganymede, and Callisto 
for the Europa Orbiter Mission, The principal goal of the (our design is to lower arrival 
V~ for the final Europa encounter while meeting all of the design constraints. Key 
constraints arise from considering the total time of the tour and the radiation dosage of a 
tour. These tours may employ J4 or more encounters with the Jovian satellites, hence 
there is an enormous number of possible sequences of these satellites to investigate. We 
develop a graphical method that greatly aids the design process. 


Introduction 

The Europa Orbiter Mission is currently scheduled to 
arrive at Jupiter by the end of the decade. The mission 
will investigate the possibility that liquid oceans may 
exist beneath the surface ice of Europa. It will attempt to 
map these regions of liquid water for follow-up missions 
to Europa. The recent discovery of life in the icc of l^ake 
Vostok, a lake deep beneath the Antarctic ice cap, lends 
impetus to Europa missions with the suggestion that life 
may be possible on Europa. 1 

In order to place the spacecraft into orbit about 
Europa, the arrival must be reduced ax much as 
possible prior to orbit insertion. This paper investigates 
the problem of lowering the arrival V w with a tour (i.e. a 
sequence of gravity assists) of the Jovian satellites, 
Europa, Ganymede, and Callisto. 

This tour is only one phase of the Europa Orbiter 
mission. After arriving at Jupiter, a maneuver will be 
performed to capture the spacecraft about Jupiter in an 
orbit that encounters Ganymede. Our tours start with 
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variations of this Ganymede encounter. After the tour 
reduces the final arrival V lT> at Europa, the endgame 
begins. The endgame is designed by the Jet Propulsion 
Laboratory (JPL) to use a combination of Europa flybys, 
small maneuvers, and 3-body effects to reduce the energy 
ol the orbit further prior to Lhe orbit insertion maneuver 
(Sec Johannesen and D’Amariu) 2 . 

Constraints for Tour Design 


Table 1 Initial conditions at Ganymede 


Launch 

Period 

v„ 

(km/s) 

Periapsis 

(Rj) 

Period 

(days) 

Early 

8.18 

9 8 

200.2 

Middle 

8.47 

9.4 

199.7 

Late 

8.14 

9.8 

191.4 


We start with a set of initial conditions at Ganymede, 
which vary depending on when the Orbiter is launched 
from Earth. JPL categorizes these conditions as “kite, 1 ’ 
middle, and early launch period. Typical initial 
conditions from each launch period are included in Table 
1. In Table K launch period ranges from 11/10 to 
1 1/25/2003, ahd results in arrivals at Jupiter from 2/28 to 
12/05/2007. Storting from initial conditions such as those 
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in Table 1, we then proceed to design the tour subject to 
various constraints. 

There arc many constraints that must be met by the 
Tour. Most important is to have as low a V* at Europa as 
possible. Based on the Hohmann transfer from 
Ganymede to Europa the lowest ballistic achievable is 
1.49 km/sec. Pcnapsis of any orbit in the tour must be 
greater than S.8 Rj (Jovian radii), to mitigate the effects of 
radiation exposure. Flyby altitude at each satellite must 
be greater than 100 km at each satellite in general, and 
must be greater than 200 km during die first flyby of any 
satellite, in order to avoid crashing into the surface due to 
navigational uncertainties. While in transit between any 
two satellites, the spacecraft must not approach within 
50,000 km of any third body (i.c. a ‘‘non -targeted” flyby) 
in order to avoid perturbing the orbit too much. Another 
design constraint is to minimize the total number of 
flybys, since each flyby may require a slight correctional 
delta-V. No close flybys are allowed when Jupiter is in 
solar conjunction. It is highly desired Lhat the tour should 
be completed while the spacecraft is within 5 AU of the 
Earth. The combination of the solar conjunction 
constraint and the 5 AU constraint limits the time of flight 
for the tour to a period that varies from roughly 280 to 
500 days, depending on whether the tour is from the late, 
middle or early launch period. Each leg of the tour must 
pass through apoapsis to allow for trajectory correction 
maneuvers. Finally, each tour must end in a resonant 
orbit with Europa. 

Tabic 2 Maximum VL for a given resonance 2 3 


Resonance 

V_ (km/s) 

3:1 

3.2 

5:2 

3.6 

2:1 

3.0 

5:3 

3.1 

4:3 

IS 

6:5 

1.2 


The endgame follows the tour. The endgame consists 
of a scries of Europa flybys combined with a maneuver at 
apojove*. The maneuvers raise penjove and lower V.„ 
while the flybys reduce the period. There is a minimum 


desired for a given final resonance acliicvcd by the 
lour 3 , as shown in Tabic 2. For example, for a 4:3 
resonance (4 spacecraft revs : 3 Europa revs) the arrival 
V,^ at Europa should not exceed 1.8 km/sec. On the other 
hand, a 6:5 resonance requires a of less than 1.2 
km/sec. which is not achievable hallistically. Since it is 
possible to achieve less than 1.8 km/sec at Eurupa for die 
4:3 resonance, most tours end with a 4 3 resonance. 

Solution Approach 


Tour 99-02 



Figure 1 Baseline tour form Europa orbitcr 


STOUR (Satellite Tour Design Program) is a software 
tool that was developed by JPL for tile Galileo mission 
tour design 4 . It has been enhanced and extended at 
Purdue to enable the automated design of gravity-assist 
tours in the Solar System as well as the satellite system of 
Jupiter 1 *. STOUR uses the patched-eonic method to 
calculate ail gravity-assist trajectories meeting specified 
requirements. 

We use STOUR as our principal tool for the design of 
Europa Orbitcr tours. From a starting condition at 
Ganymede, STOUR finds trajectories for a given path, i.c. 
a sequence of gravity-assist bodies. The massive number 
of trajectories produced by STOUR must he sifted 
through to find viable tour candidates. 

Tour 99-02 (the second tour we discovered in 1999) 
uses 15 flybys of Europa, Ganymede, and Callisio (Figure 
1)H Even with the initial position and velocity specified, 
there are tens of thousands of possible tours lhat follow' a 
specified path. The calculation of these can take weeks 
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lor a single path. When wc consider that there arc 3 n ( 1.6 
million] possible paths that begin at Ganymede and reach 
Europa in 15 flybys, we sec that the problem ol 
calculating all possible tours is intractable with current 
computer technology. Clearly, wc need to know what 
paths have the most promise to yield viable tour 
candidates before even beginning an STOUR run. 

We began tackling this problem by choosing paths by 
trial and error tempered with engineering judgement. For 
instance, we could lower the spacecraft's period and thus 
decrease the total energy relative to Jupiter in an attempt 
to reduce the final arrival V*. at Europa. A series of 
pump-downs with Ganymede would accomplish this 
quickly, but would also quickly lower the periapsis into 
the hazardous radiation environment (i.e. fry the 
spacecraft). We could then modify this path to include 
Europa and Callislo. Following such logic, we found that 
although Europa has less gravity to assist us, it is able to 
reduce period more than Ganymede for the same decrease 
in periapsis height. Wc also noticed that Callisto is very 
handy for raising periapsis. as it can do so with the lowest 
cost in increased Orbit period. If we combine these 
satellites in the right order (c.g. Ganymcdc-Callisto or 
Ganymede-F.uropa-Callisto), we could reduce period and 
periapsis at the end of a sequence of satellite flybys. The 
identification of useful path segments such as these look- 
months of experience with the problem. 
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Figure 2 P-r p plot 

To improve this trial and error method, we next 
conducted exhaustive searches tlirough all possible five- 
body path segments for the beginning of the tour. Even 
limiting the paths to five bodies left us 


computationally intensive and Lime consuming process 
which needed to be repeated for each different initial 
condition at die first Ganymede encounter. Moreover, the 
results of this endeavor were hard to interpret. A key 
question is how to characterize what will end up being a 
good tour after only five flybys. One figure uf merit is 
the at die fifth flyby, but it is difficult to draw 
comparisons between the' final VJs of path segments 
ending at different satellites. 

During the initial process we found that tracking both 
period and periapsis could often identify interesting path 
segments. Since the satellites we are working with arc in 
almost circular orbits about Jupiter, period and periapsis 
prescribe both the spacecraft’s orbit about Jupiter and the 
V„ at each satellite. 

This observation suggests the “P-rp’* plot (Figure 2). 
This is a plot of period versus perijove for orbits with less 
than 200 day periods that meet the perijove constraint ( > 
S.S Rj). The plot shows contours of constant V* for each 
satellite, assuming circular, coplanar orbits. A gravity 
assist rotates the V* vector of the spacecraft along one of 
these contours modifying the orbit about Jupiter. Where 
contours Jfom different satellites intersect, there exists a 
transfer between those satellites. These contours rive the 
values of at each satellite for this transfer arc This 
provides a method for comparing the at different 
bodies. 

If we constrain the flybys to have a minimum altitude 
of 100 kin above the Surface of the satellite, we are 
limited in how far wc can travel along a contour in one 
flyby. This is illustrated on the plot by tick marks. From 
One tick mark On a contour wc may move a maximum of 
the distance to the next tick mark up ur down that contour. 
The tick marks also can help us judge how far one flyby 
can move up or down a contour even when not starting 
from a tick mark. 

We can now see on the plot in a few minutes what 
before took months. Remembering that our goal with the 
tours is to decrease the spacecraft’s period but still keep 
the periapsis high, wc can see that Europa is most 
effective in lowering period with a minimal cost in 
periapsis height by the slope of its V*, contours. 
However, due to the distance between the tick marks, 
Ganymede is much more effective in lowering period 
with a single flyby. The slope of Callisto’s contours show 
that it is the best choice for raising periapsis as it costs the 
least in terms of increased period to do so. 


with 
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With one of these charts and a pencil, a tour designer 
can quickly sketch out a promising path lor analysis in 
STOUR. Also, known tours can be plotted and examined 
for possible improvements. 

The P-r p plots can be derived from Tisserand’s 
criterion. Tisserand showed in the 19 th cenmry UJ that 
comet orbits perturbed by Jupiter’s gravity will satisfy 
Jacobi's integral. The resulting equations can he .solved 
and plotted cm a P-r p plot. We used Tisserand 's criterion 
to verify our P-r p plots. 
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Figure 3 Hohmann transfers 

Figure 3 shows the Ganymcde-Europa, Callisto- 
Ganymede, and Callisto-Europa Hohmann transfers. 
These orbits provide a lower bound of 1.49 km/sec lor 
arrival at Europa. The chart shows this can only be 
achieved via multiple Ganymede-Europa arcs at the end 
of a tour as opposed to directly after a Cullisto flyby. 

Wc arc currently extending this method to search for 
the fastest possible path to Europa as well as low radiation 
paths to Europa. This involves developing software to 
automatically traverse the P-r p plots to find possible paths 
and calculate a cost for those paths. 

Results 

Altogether, we discovered 35 tours in 1999. Tour 99- 
02 (see Figure 1) is currently being used as a baseline by 
JPL. The details of Tour 99-02 are in Table 3. 

Tour 99-02 is one of our earlier tour designs, where we 
relied primarily on trial and error to discover and link 
promising path segments. A good example of such a 
segment is the first 5 flybys of Torn* 99-02. We start out 


with 3 Ganymede resonances, followed by a Europa- 
Callisto combination. This pattern of multiple Ganymede 
flybys followed by a Europa-CalJisto pairing accounts for 
the great majority (19) of tours w c discovered for the 
early launch period initial conditions. For low radiation 
tours, we would like the periapsis to remain as high as 
possible 


Table 3 Tour 99-02 summary 


Event #/ 
Satellite 

(km/s) 

Period 

(days) 

(Rj) 

Time 

(days) 

1 /Ganymede 

7.85 

64.3 

10.3 

0 

2/Gfinymcde 

7.85 

35.7 

9.6 

64 

3/Ganymede 

7.86 

21.4 

8.6 

100 

4/Ganymedc 

7.86 

27.8 

9.1 

122 

5/ Europa 

5.11 

20.4 

9.0 

151 

6/Call isto 

6.39 

23.) 

9.8 

169 

7/Ganymedo 

7.10 

16.7 

9.1 

193 

8/Europa 

4.74 

17.7 

9.1 

211 

9/ Europa 

4.73 

16.5 

9.1 

229 

10/Callisto 

5.75 

22.0 

11.4 

247 

1 1 /Ganymede 

5.85 

14.3 

10.3 

268 

12/Ganymcde 

5.85 

10.8 

9.3 

282 

1 3/Europa 

3.34 

10.6 

9.3 

303 

14/Europa 

3.31 

8.8 

9.1 

313 

15 /Europa 

3.29 

7.1 

8.9 

331 

1 6/Europa 

3.28 



338 


An orbit with a periapsis above 12 Rj essentially does 
not contribute to the radiation hazard 3 . The periapses in 
Tour 99-02 never exceed 12 R Jt and are rarely greater 
than 10 Rj, because when we designed Tour 99-02 
radiation was not an official constraint. The flybys of 
Europa on events 8 and 9 appreciably increase the 
radiation dosage of Tour 99-02. Since Europa has a semi- 
major axis of approximately 9.4 RJ, any flyby of Europa 
will have a significant radiation dosage. For this reason, 
our later tours avoid going to Europa until the c.nd ol the 
tour. However, the early flybys of Europa in Tour 99-02 
do serve a purpose. A glance at the P-Rp plot (Figure 2) 
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wili confirm that Europa can efficiently pump down the 
orbital penod with only a slight lowering of the periapsis. 
Tour 99-02 achieves a final V., of 3.28 krn/sec, which 
meets the maximum constraint of 3.50 km/sec imposed hy 
JPL. Later tour efforts achieve lower V„, but at a cost in 
time of night. 


Table 4 Tour 99-35 summary 


Event #/ 
Satellite 

V, 

(km/s) 

Period 

(days) 

Perijove 

(Ri) 

Time 

(days) 

1 /Ganymede 

5.99 

50.1 

12.5 

0 

2/Ganymede 

5.99 

30.5 

11.9 

50 

3/Call isto 

6.31 

41.9 

13.5 

S4 

4/Ganymede 

4.93 

21.5 

12.6 

124 

5/Ganymede 

4.93 

13.3 

11.4 

145 

6/CaUisto 

3.93 

18.0 

14.9 

155 

7/Ganymede 

2.37 

10.7 

13.9 

194 

S/Ganymedc 

2.37 

7.2 

11.7 

215 

9/Ganymede 

2.37 

5.5 

9.1 

222 

10/Europa 

2.45 

5.2 

9.0 

232 

1 1 /Ganymede 

1.59 

5.3 

94 

245 

l'2/Europa 

1.64 

4.7 

9.33 

253 

I3/Europa 

1.62 
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We used P-r p plots to discover Tour 99-35. First, a 
promising path lor the tour was selected from the P-Rp 
plot and evaluated interactively (in STOUR) to test its 
effectiveness. We used this run in conjunction with Lhe P- 
Rp plot to adjust our selected path as necessary. Finally, 
the selected path was used as Lhe basis of an automated 
search in STOUR. Tour 99-35 is listed in Table 4. 

WiLh Tour 99-35, we were trying to limit the number 
of flybys and maintain a high periapsi? for low radiation. 
Consequently, we started Tour 99-35 with the highest 
periapsis possible. This turned out to be a periapsis of 
13.2 Rj for an inilial condition from the late launch 
period. The use of the P-r P plot paid off nicely, as Tour 
99-35 has the lowest time of flight of any tour we 
discovered. Tour 99-35 also has a final arrival V w of 1 .62 
km/sec, which is about as close to the Hohmann limit of 
1.49 km/sec as we have been able to achieve. The 
radiation environment during Tour 99-35 is excellent 
through event 10, and if we had stayed at Europa on event 


10. Tour 99-35 would have an exceptionally low radiation 
dosage. However, we chose to tack on an additional 
Ganymede-Europa sequence to lower the final arrival V. A 
from 2.45 km/scc to 1.62 krn/sec (a considerable 
improvement). Consequently, we take a hit in radiation 
dosage on events 1 1 and 12. 

The respective paths of Tours 99-02 and 99-35 appear 
in Lhe P- r p in Figure 4. A comparison of these tours 
demonstrates the efficacy of the P-r p plot. The path of 
Tour 99-35 is represented by dashed lines; that of 99-02 
by solid lines. From the point of view of path selection, 
we can see a clear inefficiency in Tour 99-02 for events 
G3 and G4. G3 pumps all the way down to a periapsis of 
8.7 (which is a slight violation of the r p constraint), and 
then G4 pumps up to a transfer to Europa (E5). Instead of 
Lhis roundabout method of reaching E5, in retrospect we 
could have simply used the G3 transfer to reach Europa, 
thus saving a flyby and reducing the radiation dose. A 
similar inefficiency for Tour 99-02 occurs with the E8 


ip* 


Tour 99-02 


Ci 


Cl 


Tour 99 35 


G2 


G2 


<? g4 

-sr 

es y *' 

^^ClO — - 

E12 } r GH 
E13 


r:<t 


EX0 

Gt 1 


Gif 


G9 


G7 




Figure 4 P-r p comparison of Tours 99-02 and 99-35 


and E9 flybys. On the other hand. Tour 99-3 5 proceeds 
smoothly from initial condition to final arrival. There is 
very little "wasted movement” or meandering about the P- 
r p plot. Furthermore, in general each flyby in Tour 99-35 
moves farther along a V w curve than the flybys of Tour 
99-02, implying more efficient use of each flyby. Thus, 
the use of the P-r p plot helps us in two particular ways. 
First, we can always select the next best event in a path. 
Second, since we have a good way of selecting a path, 
each flyby can be more effective (move farther along a 
curve). In addition, the P-r p plot can also be used to 
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Figure 5 P-r p plot of initial conditions 


evaluate the efficacy of an existing tour at a glance. The 
result, as can clearly be seen in Figure 4, [$ a shorter, 
more efficient tour. 

Tour 99-35 also benefits from having a better initial 
condition. When we designed Tour 99-35, our goal was a 
lower radiation dosage, so we selected the highest initial 
rp that we had available. All the initial conditions that wo 
used for tour design are plotted in the P-r p plot in Figure 
5. The beginning launch period initial conditions are 
marked with n ‘TP, the middle launch period initial 
conditions with an “M* , and the late launch period initial 
conditions with an **L'X By comparing the positions on 
the P-r P plot of the beginning, middle, and late launch 
periods, we can see that the initial GGEC sequence for 
most of our beginning launch period tours intuitively 
makes sense. Given the low periapsis, we need to pump 
down quickly to gel to one of the flatter Cailisto V* 
curves, which are efficient at increasing periapsis while 
not increasing period too much. We can also see that 
since many of the middle and into launch periods start 
with higher r p values, we can reduce period somewhat 
more at Lhe beginning of the tour without lowering the 
periapsis too much (and thus our radiation dosage is much 


lower). (Of course, there is a delta- V cost associated with 
starting the lour at a lower periapsis). Also, the time of 
flight for the late and middle launch periods is generally 
lower, because they start with a lower V* value. Clearly, 
the initial conditions greatly affect our tour design 
Strategy. 

Tabic 5 lists the best tours for arrival time of 
fiight, radiation dose, find number of flybys for each 
launch period. The table does not include tours that had 
better performance buL violated either the flyby altitude, 
non-targelcd, or solar conjunction constraints One issue 
that has not been addressed in detail is phasing (ie, timing 
between the jovian satellites and the spacecraft orbit) 
Phasing can significantly affect the performance of a tour 
The P-r p plot does not address phasing or liming issues at 
all. Hence, while a path may look promising on the P-r p 
path, STOUR may not be able to find viable lours 
Phasing tends to become more of an issue towards the end 
of the tour, when the energy of t lie spacecraft orbit is low 
and fewer transfer arcs between satellites are possible. 

Given the combination of V w and low radiation 
constraints, we almost always want our last Callisio- 
Ganymede transfer orbit to have a periapsis as close to 
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Ganymede's semi-major axis (14.971 Rj) as possible 
(since we arc trying to achieve a Hohmann transfer 
between CaiJisto and Ganymede). In practice, due to 


Table 5 Details of Best Tours 


Parameter 

Early 

Middle 

Late 

1 .Arrival Date 

99-02 

99-23 

99-35 


(2/01/08) 

(7/27/OS) 

(S/1 3/OR) 

2. Lowest 

99-33 

99-25 

99-35 


0.59) 

(1.71) 

(1.62) 

3. Lowest Dose 

99-18 

99-26 

99-35 


(9.2) 

(S3) 

(7.2) 

4. Fewest Flybys 

99-04 

99-32 

99-35 


(14) 

(13) 

(12) 


phasing, the ideal transfer between Ganymede and 
Callisto proves elusive, as docs the later final Ganymede- 
Furnpa transfer. In fact, our experience is that the final 
sequence of flybys is much more of a limiting factor than 
any other portion of the tour (ie. in the middle of the tour, 
many transfer orbits for a given flyby are available, at the 
eud, only a few). 

Conclusions 

The process of designing tours for the Europa Orbitcr 
mission has been considerably streamlined via the use of 
die P-r p plots. Delta V savings exceeded expectations. 
The V. improved from 3.5 km/sec for the initial tour to 
1.6 km/sec lor the best case. Radiation dosage lor the 
later lours was reduced up to 70% from the earlier 
numbers, and exceeded the expectations of the radiation 
group at JPL n . P- r p plots help tremendously with path 
selection and insight, and have proven a useful tool in 
other gravity assist missions. In the process of designing 
and improving the Europa Obiter mission, we found 
something with a general applicability. 
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